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ABSTRACT 

It  has  been  observed,  that  the  behavior  of  missiles  either 
■ovine  under  the  Influence  of  non-linear  forces  and  *a*entB  or  flying 
at  large  angles  of  yaw  Is  frequently  veil  described  by  curves  of  the 
earns  for m  as  those  generated  by  linear  force  systems  and  small  a-Kles 
of  yaw.  With  this  In  Mind  an  "equivalent"  Unear  solution  to  the 
actual  equations  of  yawing  motion  is  obtained. 

This  equivalent  linear  solution  has  been  used  in  the  analysis 
of  a  vide  variety  of  program  fired  on  BRL's  Spark.  Ranges  and  considerable 
success  ha-,  been  experienced.  Excellent  internal  consistency  has 
been  observed  in  Measuring  non-linear  nomal  and  Magnus  forces  and 
their  moments  and,  in  all  cases  where  vind  tunnel  results  vere  available, 
they  vere  in  good  agreement  with  range  results. 

The  application  of  this  technique  to  the  equally  important  problem 
of  predicting  yawing  Motion  is  described. 
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(u^,  Ug,  Uj)  velocity  vector 
1*1  magnitude  of  the  velocity 
magnitude  of  the  velocity  at  mid-range 
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*10  -  9  ■  *  -  *20 

arccos  n^ 

(°r  Qg,  Qj)  angular  velocity  vector  of  co-ordinate  system 
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circumflex  superscript  with  the  exception  of  v  denotes 
quantities  In  the  fixed-plane  co-ordinate  system 

2 

tilde  superscript  ^rith  the  exception  of  8  denotes  quantities 
appearing  in  the  solution  of  the  li  nearized  yaw  equation 

asterlak  denotes  the  modification  oi  "usntitles  involving 
aerodynamic  coefficients  through  the  consideration  of  the 
cosine  of  the  yaw  angle,  t. 


1. 


INTRODUCTION 


The  free  flight  spark,  range  technique  measures  the  aerodynamic 
forces  and  moments  acting  on  a  missile  by  means  of  very  accurate 
observations  of  its  motion  in  flight .  This  process  requires  a 
knowledge  of  the  functional  dependence  of  these  forces  and  moments 
on  the  dynamic  variables  of  the  motion  in  order  that  the  solution 
curves  to  the  differential  equations  of  motion  may  he  obtained. 

These  solution  curve«  are  fitted  to  the  motion  and  the  forces 
and  moments  calculated  from  the  parameters  of  the  fit. 

This  need  for  solutions  in  closed  form  has  traditionally 
limited  range  tests  to  motions  which  are  described  by  linear  equations. 
Since  non-linear  terms  arise  from  both  the  size  of  the  motion  (non¬ 
linear  geometry  terms)  and  the  presence  of  second  order  or  higher 
terms  in  the  aerodynamic  force  expansion  (non-linear  force  terms), 
this  raeans  that  the  range  technique  is  restricted  to  configurations 
possessing  linear  force  systems  and  flying  at  small  angles  of  yaw. 

Strangely  enough,  a  number  or  models  have  been  fired  in  the  BHL 
R antes  which  either  possessed  known  non- linearities  In  their  force 
systems  or  flew  at  large  angles  of  yaw  and  it  was  found  that  their 
motion  could  be  very  well  fitted  by  functions  which  were  solutions 
of  the  linearized  equations.  This  seemed  to  imply  that  the  parameters 
of  these  linear  equations  should  be  "average  values"  of  the 
coefficients  of  the  parent  non-linear  equation.  It  is  the  purpose 
of  this  report  first,  to  derive  the  equations  which  relate  these 
average  values  to  the  non-linear  force  terms  and  certain  characteristics 
of  the  motion  and  then  to  demonstrate  the  great  value  of  these 
relations  by  applying  them  to  s  number  of  programs  which  have  been 
fired  at  BRL.  The  success  of  this  technique  more  than  doubles  the 
values  of  free  flight  ranges  for  both  the  ballistlcian  and  the 
aerodynamicist . 

Finally  the  extenaior  ’.his  method  to  the  even  more  Important 
problem  of  the  prediction  01  j  ing  motion  is  described. 
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i.  ikwjgc  cou'imx/noua 

In  order  t  o  investigate  the  yawing  motion  of  miss  ties  moving  at 
large  angle:;  of  yaw,  the  c/.u:/t  equations  of  motion  raur-t  be  derived. 

An  Important  feature  of  this  derivation  lies  in  she  proper  selection 
of  the  co-ordinate  system  and  the  dynamic  variaries .  This  selection 
should  be  if"- do  go  that,  the  resui tirg  equations  ore  no  simple  as 
possible  and  '•easonnbly  compatible  with  the  basic  assumption  tliat  their 
bolutior.  may  be  approximated  by  a  solution  of  the  linearised  equations. 

In  other  words  the  non-linear  equation  must  possess  the  same  general 
characteristics  of  the  linearized  equation,  i.e  ,  order  of  enuatlors, 
number  of  variables,  symmetry  of  variables,  etc. 

For  arbitrary  rigid  bodies  the  best  co-ordinate  system  is  one 
whose  axes  lie  along  the  body’s  principal  axes  of  Inertia.  Since  toe 
mass  distributions  of  cost  missiles  are  rotatlonully  symmetrical  about 
the  longitudinal  axis  of  inertia,  all  transverse  momenta  of  inertia 
are  assumed  to  be  equal  and  all  axes  perpendicular  to  the  longitudinal 
axis  are  principal  axes  of  mertid,  A  right-handed  Cartesian  co-ordinate 
sy8tem  with  numbered  axes  will  be  constrained  so  that  its  l-axis  is 

aligned  with  the  missile’s  longitudinal  axis  and  pointing  forward. 

■* 

fc  ,  the  angular  velocity  vector  of  the  co-ordinate  system 
relative  to  on  inertial  system,  will  have  components  (fl  ,  fl-,,  fl,) 

X  c.  J 

while  the  missile’s  angular  velocity  vector,  to,  will  have  components 
(a^,  (h,,  aij).  This  definition  of  the  co-ordinate  system,  then,  requires 
that  «  uv,  and  flj  «  a>y  If  and  an  initial  orientation  of  the  2- 
axls  are  specified,  the  co-ordin.ate  system  would  be  completely  determined. 
In  Appendix  A  the  equations  of  motion  are  derived  for  arbitrary  values 
of  Oj  and  we  see  that  these  equations  would  be  greatly  simplified  if 
«  0.  For  this  reason  we  will  maJte  considerable  use  of  a  non-spinning 
co-ordinate  system  for  which  the  2-axl3  initially  lies  in  the  horizontal 
plane  pointing  to  the  left  and  whose  axial  spin,  fl^,  in  identically  zero. 

Turning  now  to  the  question  of  dependent  variables  it  may  be  seen 
that  the  major  contenders  are  the  Eulerian  angles  and  t»-e  direction 
cosines.  The  question  of  ccmpatability  now  appears •  Althougu  toe 
linearized  equations  in  terms  of  Eulerian  angles  are  symmetric  in  the 


components  of  yaw,  the  corresponding  exact  equation*  cannot  be. 

This  is  clear  fiwj:  chc  definition  of  the  angles.  For  this  reason* 
we  will  express  c-_-  equations  in  terms  of  direction  cosines. 

At  this  point  ve  must  precisely  define  our  variables.  The  yaw 
angle  is  defined  to  be  the  angle  between  the  missile’s  axis  and  the 
tangent  to  its  trajectory.  If  the  missile's  velocity  vector,  u, 
has  components  (u, ,  Up,  u^)  and  magnitude,  u,  the  sine  of  the  yaw- 


angle  is’-— 


and  a  unit  vector  lying  in  the  plane  of  yaw  has 


direction  cosines  (0, 


J  2  2 

^  +  u3  ^“2  +  u3 

consider  only  those  missiles  possessing  trigonal  or  greater  rotational 
symmetry,  it  will  he  convenient  to  represent  quantities  in  the  plane 
normal  to  the  missile's  axis  by  complex  numbers.  With  this  in  mind 
we  define  the  complex  yaw  vector,  X.,  to  be  a  vector  which  lies  in  the 
plane  of  yaw  and  whose  magnitude  is  the  sine**of  the  yaw  angle. 


-).  Since  we  will 


The  cosira  of  the  yaw  angle  will  he  needed  in  this  report  and  will  he 
denoted  by  i. 


As  can  he  seen  from  an  examination  of  Reference  1  the  julerian 
eagles  also  introduce  considerable  algebraic  complexity. 

» 

In  Reference  2,  is  used  as  e  characteristic  velocity  for  forming 

dimensionless  quantities  and,  hence,  the  cceplex  yaw  of  that  re¬ 
ference  is  proportional  to  the  taagert  of  the  yaw  angle.  Although 
this  use  of  does  simplify  the  center  of  miss  relations,  it  else 

introduces  a  number  of  complications  which  oui  choice  of  u  avoids. 


%’ikV  \iirri4* 


In  a  similar  fashion  the  axial  ana  transverse  component;  of  the  ang  lar 
veloclty  nay  be  separated  and  written  in  complex  form. 


(<xu  +  l<o,)d  (3) 

U  m  .  - .... - 

U 

where  A  is  diameter. 

Finally  we  will  select  our  independent  >crlablr  to  be  distance  measured 
along  the  trajectory  in  calibers.  If  this  variable  is  identified  by  p. 


where  t  is  tine. 

The  generalized3  lines'-  er.wiiUi  jf  the  aerodynamic  force  and 
went  assumes  that  tb*  f-rce  and  women*,  are  linear  functions  of  X,  u, 
and  the*-.  d-;i  .-.tiyes  in  a  non-rolling  ;o- ordinate  nystea.  For  syMaetric  / 
•..b_.^.es  this  assumption  introduces  ei  jhteen  coefficients.  If  we  Unit 
ourselTSs  to  onTy  those  coefficients  isring  a  measurable  effect  on  the 
Motion,  this  total  reduces  to  ten  ar-.  we  have  the  following  expansion.  / 

Fx  -  -  p  d2  u2  (?) 


f2  +  if3 


-  pd3u2  vK 


d2  u2  £  (  Kjj  +  iv!Cy)X  +  iKgU  -  K^iX  +  it>  •  ]  J 

vKA  <7)/ 

3u2[(  -vl^  -  1Kjj)X  -  KgM  -lK^d't  HxA 


*2  ♦  i*j  -  Pd3u2|_(  viqj  -  lK^X  -  Kf  -iK^(x\  i*X)  : 

i  (F^,  »2,  Fj)  are  -  cnpooents  of  the  aerodynaadc  fore  L  > 
(Kj,  Ug,  Kj)  are  components  of  the  aerodynamic  nor/^t, 
p  is  air  danalt-  . 

¥  .  ,  end 

f 

r  are  Alans/ lodees  aeroAnm?  coefficient*,  r 


'--fr ■ 


For  this  report  ve  will  only  consider  non-linearities  caused 

by  the  size  of  the  yaw  angle.  Since  symmetry**  requires  that  the 

—  ■  2 

aerodynamic  coefficients  must  be  functions  of  XX  =  &  ,  we  will  only 
consider  such  a  dependency  on  the  square  of  the  sine  of  the  yaw  angle. 

Two  other  force  coefficients  will  be  used  in  this  report.  These 
are  based  on  a  resolution  of  the  aerodynamic  force  in  the  plane  of  yaw 
:  ’.to  components  along  the  trajectory  and  perpendicular  to  it.  If  FQ 

ii  the  component  along  the  trajectory  and  FT  the  component  perpendicular 

X* 

to  the  trajectory  and  pointing  toward  the  missile's  nose. 

FD  -  -  Pd2u2  Kjj  (9) 

Fl  -  -  pd2u2  S  (10' 

-here  8  is  the  magnitude  of  the  sine  of  the  yaw  angle.) 


In  order  to  obtain  relations  between  these  two  new  coefficients 
and  those  defined  in  Equations  (5  -  6),  we  require  that  the  2-axis 
be  in  the  plane  of  yaw  so  that  X  is  real  and  equal  to 


JFi-8fa 


(ii) 


?L  -  8  *X  +  *  *2 


(12) 


I 


Substituting  Equations  (5),  (6),  (9),  and  (10)  in  (11)  (12) 

and  neglecting  the  Kg  and  terms. 
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m  is  mass. 


J.  SMALL  YAWING  MOTION  WITH  NON-LINEAR  MOMENT  AND  NO 

CAMPING 

In  this  section  a  verv  simple  example  of  non-linear  yeving  motion 
will  he  considered  in  some  detail .  This  is  the  case  of  a  missile 
flying  at  small  yaw  and  acted  on  by  a  cubic  overturning  cr  restoring 
moment.  Later  the  effect  of  other  non-linear  forces  and  moments  as 
well  as  the  geometric  non-linearities  will  be  considered. 

The  basic  feature  of  all  the  theoretical  results  of  this  report 
is  tne  assumption  that  over  "small"  sections  of  a  missile's  trajectory 
the  non-linearities  present  in  the  equations  of  motion  do  not  cause 
the  motion  to  be  qualitatively  different  from  motion  based  on  the 
linearized  equations.  For  example,  the  yawing  motion  of  a  symmetric 
missile  acted  on  by  non-linear  forces  and  moments  should  still  bo 
epicyclie  when  a  “small"  portion  of  the  trajectory  is  considered.  The 
parameters  of  the  epicycle,  however,  would  pr.Lably  be  related  to  the 
size  of  the  motion.  This  assumption  seems  to  be  reasonable  when  the 
non-linearities  are  themselves  “small”.  The  experience  of 
balllstlcians  as  based  on  actual  free  flight  tests  have  indicated  that 
relatively  long  section;;  of  trajectories  and  large  non-linearities  are 
still  “small"  enough  for  this  assumption. 

In  order  to  illustrate  this  point  the  data  analysis  of  a  non¬ 
linear  crag  force  will  be  outlined.  Since  this  is  the  only  non- lineal- 
aerodynamic  force  which  un  to  the  present  time  has  been  .accessful ly 
handled  by  ballicticians,  this  digression  will  also  provide  a  good 
background  for  later  derivations. 

The  non-linear  dependence  of  the  drag  force  on  the  magnit’ide  of 
the  yaw  is  very  well  described  by  the  assumption  that  the  drag  coefficient 
is  a  quadratic  function  of  8  where  8  is  the  aim  of  the  total  yaw  angle. 


I  e  • 


In  Reference  5  it  is  shown  that  the  usual  drag  reduction  for  flat 
trajectories*  involves  the  fitting  of  the  time  measurements  to  a 
cubic  in  distance. 

2  x 

t’‘aD  +  alp+B2p  +a5p  (22) 

where  p  is  measured  from  the  middle  of  the  observed  trajectory. 

It  is  further  shown  there  that  under  certain  reasonable  approximations 
the  following  equation  applies  for  drag  force  of  the  same  form  as 
Equation  (21). 

<25) 

where  is  velocity  at  midrange  (p  =  0), 

is  a  constant  related  to  drag  dependence  on  Mach  number, 

t  ?  2 

I*  (p)  *  V  \  &  dr  d<k  and 
o  o 

J_  -  K_ 

Di  *  \ 

If  the  yaw  is  veil  described  by  an  epicycle,  then 

<*> 

-  aj 

where  Kj  -  e  J  -  0^  +  0j'P»  aaA 

KJ0*  ay  ^JO*  ^JO  ”*  real  coaBtants. 

2—22  “  flo)  i(0g  “  0j) 

8  «  XX  •  Kj_  ♦  ^  (e  A  +e  )  (25) 

Although  the  expression  in  parentheses  may  be  given  more  simply  by 
2  cos(01 02),  this  form  will  be  more  convenient  for  the  non-linear 
yaw  equation.  For  the  case  of  zero  yaw  drag  coefficient,  we  see,  from 
Equations  (22)  and  (23),  that 

s  ^  (a6) 


A.  trajectory  is  said  to  be  flat  when  the  component  of  the  gravational 
force  along  the  trajectory  does  not  change..  The  yaw  drag  treatment  may 
be  extended  to  treat,  non-flat  trajectories  If  this  is  necessary. 


When  the  yew-drag  coefficient  has  a  censurable  effect,  it  does  not 
usually  affect  the  quality  of  this  cubic  fit  for  segments  of  trajectories 
as  long  as  10,000  calibers.  It  has  been  found  that  in  this  case  tne 
same  configuration  flying  at  different  average  yaws  will  provide  good 

^a2 

cubic  fits  but  dlffere  si  value:  of  -  .  The  problem  is  then  to  find 

B1 

the  quadratic  contribution  of  the  yaw-drag  term  of  Equation  (23). 


Since  all  drag  reductions  contain  a  number  of  periods  of  the 
cosine  term,  the  cubic  fit  of  the  time  history  effectively  neglects 
this  component  of  &  .  With  this  in  mind  we  look  for  the  "best"  quadratic 
approximation  to  I, (p)  over  the  length  of  observed  trajectory,  L. 

r~> 

2 

To  do  this  we  define  an  average  squared  yaw,  6  ,  so  that  the  Integral 


ft 


^el  T?  2 

(p)  -  8  (•£-)  dp  is  a  minimum.*  6  is,  therefore,  the  best 


quadratic  coefficient  from  the  standpoint  of  least  squares.  Differentiating 
the  integral  with  respect  to  8  and  setting  the  result  equal  to  zero. 


mj  I* 

.2  -5 


L 

^  p2  Ix(p)  dp 


2 


(l)p4  dp 


that  this  average  squared  yaw  is  different  from  the  more  tradition¬ 


ally  used  mean  squared  yaw  ? (K^2  +  Kg2)  .Ip.  Although  the  average 


squared  yaw  has  better  theoretical  Justification,  the  numerical  difference 
is  usually  not  important. 
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where  I^p) 


-  SOgP 


„  2|  *  ‘  1  +  ^  I  .  „  2  e  -  1  +  20gP 

'H  J  “  L  «?■  _ 

2 

The  complete  expansion  for  8  may  be  computed  from  the  power  aeries 
expansion  of  I^(p) •  Since  three  non-i-aro  terms  are  usually  sufficient, 
the  following  expression  for  8“  is  easily  derivable. 

"2  r  2 r  1  (aiL)2  1 

8  “  10b°  |_5Pr>  +Tnrn  +^rsT)  J  (27 ) 

2  x  (o^hf  (a^)4 1 1 

+  *ao  |*2T)  +  THTH  +1T5T)J  j 

Returning  to  equation  (26), 

h  -  V®2  (28) 

range  pd  T.  8  ^ 

2 

In  most  cases  a  good  approximation  fcr  8  can  be  obtained  by  taking 
the  first  term  in  the  expansion  of  Eq.  (27*). 

•*•*!>  “  %>  +  V  (K1°2  +  ^02)  (29/ 

range  o 

The  value  of  this  relation  is  clearly  shown  in  Figure  1  where 
Kjj  is  plotted  versus  8^  for  a  body  of  revolution  fired  at  yaw 
•ngSPup  to  50°.  A  further  check  on  this  technique  lies  in  the 
good  agreement  of  the  wind  tunnel  value  of  X^2  *  1-85  with  the  free 

flight  value  of  1.82.  (See  Table  IV  on  Page  45. ) 


*  1.85  with  the  free 


Free  flight  tests  have  shown  that  the  situation  for  the  yawing 
notion  is  quite  similar  to  this  yaw  drag  case.  Although  the  non- 
linearities  do  not  measureably  change  the  nature  of  the  epicyclic 
motion  which  is  predicted  by  the  linear  theory,  they  do  affect  the 
values  of  the  epicyclic  parameters.  It  is,  therefore,  necessary 
to  obtain  relations  of  the  sane  type  as  Eq.  (29). 

Since  the  simple  case  of  the  small  amplitude  /swing  notion  of  a 
missile  acted  on  by  a  cubic  static  moment  will  first  be  considered. 
Equation  (20)  can  be  considerably  simplified.  The  3mall  amplitude 
assumption  implies  the  approximations  1=1  and  1*  *  0  while  the 
restriction  to  a  static  moment  eliminates  the  terms  in  H,  J£  and  T. 
The  equation  under  consideration,  therefore,  becomes 


X"  -  ivX*  -  (Mq  +  \  o  )X 
where  ^  ^ 


"2  m£r\z 

KHm  ^  0  + 


62  -  XX. 

The  solution  to  the  linearized  fora  of  Equation  (30)  for  a 
gyroscopically  stable  missile*  is  an  epicycle  without  damping. 

£ 

\  _  V  .  A  j  V  «  6 


*  A  missile  1b  gyroscopically  stable  if  v  >4*^ .  All  statically  stable 

missiles  (Mq<0),  for  example,  are  gyroscopically  stable.  Statically 
unstable  missiles.  (It  >0)  are  gyrosooplcally  stable  when  their 


r(l^>0}  are  gyroscopically  stable  when  their 


gyroscopic  stability  factors,  s  v  ,  are  greater  than  unity. 

TJT 

0 
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where  *  '  0  =  01O  +  f't  P 

0X  P2> 

*10’  ja^Q,  and  0^  are  real  constants. 

In  Equation  (31)  the  K, A  and?'.-  s  depend  on  the  initial  conditions  hut 
the  frequencies,  0,  arc  functions  fo  the  coefficients  of  the  linearized 
equation  and  are  Independent  of  the  initial  conditions.  The  veil  known 
relations  for  the  frequencies  are 

7X  *  *  (32) 

71  •%  -  %  (33) 


It  will  be  shown,  however,  that  the  frequencies  of  the  quasi-linear 
solution  of  Equation  (30)  do  defend  on  the  Initial  conditions. 

If  the  solution  to  Equation  (30)  is  assumed  to  be  of  the  same 
fora  as  Equation  (31),  it  can  be  substituted  in  Equation  (30)  to 
provide  the  following  equation 

*10  •  1  *  *'l  ♦  *  <  '  *0  '  **2  ^LO2  +2K202) 

♦  *20  ^  2  -  +  *  ‘  ^  <*20  +  ^lO2)  <»> 

T  a  ^i-^a)  ?  ^2-^x1 

'^[*10  *20«  +  *10*20  «  j  -°- 

*  The  tilde  superscripts  are  used  to  eaphasics  that  the  7/ a  appear 
in  the  solution  to  the  linearized  equation. 
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Since  we  know  that  the  epicycle  is  a  good  description  of  the  actual 
motion,  the  third  term  in  Equation  (3*0  which  contains  mixed  frequencies 
is  neglected  and  the  following  two  equations  may  be  obtained. 


*1  •  t  ♦  %  4 

*2  <?  -  \  4 


vh«r«  5  _  ■ 


2  2 
0  +  ^0 


4  ■  v  * 

Unfortunately  Equations  (35  -  36)  are  both  quadratic  equations  and  two 
different  values  of  each  frequency  are  possible.  This  difficulty, 
however,  can  be  resolved.  Since  we  are  considering  solutions  close 
to  the  solution  of  the  linearized  equation,  the  solution  of  Equation 
(35)  close  to^  should  be  selected  and  similarly  for  Equation  (36). 

This  means  that  the  larger  root  of  Equation  (35)  and  the  smaller  root 
of  Equation  (36)  should  be  used. 

o 

Equations  (35  -  36)  show  that  8el,  the  effective  value  of  the 
squared  yaw  for  the  1-th  frequency,  is  twice  as  sensitive  to  the 
amplitude  of  other  frequency  as  it  is  to  its  own  amplitude .  The  effect 
cases  from  the  cosine  term  vhich  was  omitted  from  the  yaw  drag  analysis. 

O 

for  the  yawing  motion  this  periodic  part  of  8  can  not  be  neglected. 

In  order  to  obtain  relations  similar  to  Equations  (32-  -  33), 
la  first  eliminated  between  Equations  (35  -  36)  and  then  v  is  eliminated 
between  them 


'*20 


01  .  02  -  Mq  ♦  ^  8/ 


K  4 


^2  ®el 


*L-h 


21 


*  +  01  V  -  h  *20 


At  first  glance  this  derivation  of  our  "equivalent”  solution 

seems  to  be  surrounded  by  an  atmosphere  of  expediency  and  to  be 

resting  on  empirical  assumptions.  Actually  this  is  not  true.  As  we 

shall  now  show  our  technique  hau  a  good  theoretical  background  and 

requires  a  minimum  algebraic  work  in  comparison  with  other  methods. 

PinaUjr  it  will  be  shown  that  the  technique  is  not  restricted  to  non- 

2 

linearities  proportional  to  6  but  can  be  easily  extended  to  polynomial 
2 

functions  of  6  . 

The  derivation  of  Equations  (35  -  36)  may  actually  be  considered 
the  first  step  of  an  Iteration  procedure.  It  ie  certainly  reasonable 
to  as 3 use  for  small  non-linearities  that  the  first  step  of  the  iteration 
has  the  same  form  as  the  linear  solution.  The  error  term, 

^  |V  *20  *  '  0fi)  'ho**?*  ' 

introduces  mixed  frequencies  which  are  a  characteristic  of  a  non-linear 
equation.  The  next  step  in  the  iteration  would  he  to  assume  a  solution 
at  the  form 

i02  1(20,  -  0o>  K202  -  0,) 

+  X20e  +K30e  +KU0  e 


Substitution  of  Equation  (39)  in  Equation  (30)  would  provide  four 

f 

complicated  equations  in  terms  of  the  four.  K,A  s  and  two  frequencies 

j  ’  10 

b,  and  our  error  term  would  then  be  of  the  form 

1(30,  -  20  )  i(30  -  20. )  100  -  302)  100,-50.) 

«2^e  +Age  2  1  +Aje  1  2  ♦  A^  e  *  ^  ). 

where  the  A^  are  fifth  order  combinations  of  the  K.  ^|s.  This  process 
may  then  be  further  iterated  to  yield  a  series  exptnsion  of  the  almost 
periodic  solution  of  Equation  (3C).  Fortunately  as  »e  shall  see,  the 
experimental  results  of  this  report  show  that  only  the  first  Btep  of 
this  process  is  needed. 


*  i  «> 


Another  method  of  treating  Equation  (30)  which  may  seer,  to  he  more 
elegant  than  the  method  of  this  report  is  that  of  Kryloff  and  Bogoliubof f . * 
It  will  he  shown  that  this  method  provides  the  same  results  as  our  direct 
substitution  approach  and  requires  more  algebraic  work.  Kryloff  and 
Bogolluboff  move  the  non-linear  terms  to  the  right-hand  side  of  the 
aquation  and  place  the  solution  of  the  linearized  equation  ?n  these 
terms  so  that  they  are  functions  of  the  independent  variable  p.  The 
method  of  variation  of  parameters  is  then  used  to  solve  the  resulting 
inhomogeneous  linear  equation.  In  ord-r  to  solve  the  differential 
equations  for  the  parametric  functions,  the  terms  arising  from  the 
inhomogeneous  term  of  the  original  equation  are  averaged  over  the  two 
periods  of  the  motion. 


As  an  illustration  of  this  method  we  rewrite  Equation  (30)  as 


X"  -  ivA  -  MoA  =•  Mp  62  X  =  Mg  f(p) 
where  f(p)  -  (K^2  +  SK^2)^  e~  1  +  (K^2  -  SK^2)^  e 
_  2.  P  -  %} 


Vac 


+  *20  *10  e 


The  parameters  of  this  solution  are  the  magnitudes  K^q  and  the  posse 
angles  8.  Differentiating  Equation  (31)  're  have 

»  »y»  $2  ,  »  .  %  ,  *  .  ^2 

*  “  **!  ho  •  +  %  ^0  •  +  <*k>  +  i0io  ho)e  + 

If  the  lust  two  terms  are  set  equal  to  zero, 

(K10  +  *  <*S»  +  (*«) 


*  In  Reference  7  it  is  shown  that  this  method  is  fundamentally  the  same 
as  that  of  Van  der  Pol.  They  differ  in  the  fact  that  Kryloff  and 
Bogolluboff  express  the  linear  solution  in  polar  co-ordinate  form 
and  Van  der  Pol  expresses  it  in  Cartesian  foim. 
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then  Equation  (1*1 )  iwy  be  differentiated  a£tiln  to  yield 

X"  -  -  0^  2  K1q  *  ^  -?2  2  X20  *  2  +  (K10  +  ^10  K10)? 


*  $2  (*20  *  $20  *20^ 

Substituting  Equations  (31,  U,  43)  for  X,  x\  X"  In  Equation  (40) 
and  solving  for  K^Q  +  i^0  X'in  by  use  of  Equation  (42), 

-  10L 

*  ^  *2  *  (UX\ 

ho  +  ^io  *io  -  i^7%r~  • 

If  the  right  aide  of  Equation  (44)  is  now  averag'd  over  a  period  of 

d  -  d  .  and  the  result  divided  into  real  and  1* aginary  parts, 
r2  rl' 


ho  ■ 0  r  -i 

„  ho  *  V 

*10  ■  :l  -X- 


By  By»*try, 


«*•«  . 
f  .  IV  ♦  V 

J 

Equation*  (45  -  47)  show  that  the  first  appraxiaatloaa  to  the 
solution  of  Equation  (30)  is  an  epicycle  with  no  damping  but  with 
frequencies  which  differ  frou  those  of  the  solution  of  the  linearised 
equation.  The  corrections  to  the  frequencies  are  the  7^0  •  Kiven 
Equations  (46)  sad  (46). 


+  £K 


10 


m 


In  order  to  compare  these  results  with  Equations  (37  -  38),  Equation 
(48)  and  (4-9)  are  first  added  and  then  multiplied.  After  malting 
use  of  Equations  (32)  and  (33),  the  following  relations  may  he 
obtained. 


< 


^>+M2 


2  2 
K10  -  *20 
v  r^A  r*ft 

r  °1  -% 


10 


*20 


)  (50) 

ho  ?L  ■  *20  % 

*  %  ■% 


(5U, 


^.o2  *  aW!'<lW!  ’  =*]/> 

2 - 


But  if  second  order  quantities  are  neglected.  Equations  (50  -  51)  are 
the  sane  as  Equations  (37  -  3u).  Thus  for  twice  the  work  ve  get  the 
sane  result. 


As  a  further  basis  for  our  quasi-linearization  technique  a  special 
for*  of  Equation  (30)  will  be  considered,  for  which  an  exact  solution 

Is  known.  This  case  Is  that  of  a  statically  stable  non- spinning  missile 

10 

in  planar  yawing  motion.  If  X  ie  replaced  by  6  e  where  ©  Is  the 
orientation  angle  of  the  plane  of  yaw.  Equation  (30)  assumes  the  fur* 


-  (My  +  Mg  82)B  -  0 


(52) 


where  My<0.  For  planer  yawing  notion  &  must  go  through  zero  and,  hence, 
the  two  amplitudes  must  be  equal  (K^  -  KgQ  *  K).  Finally  according  to 
Equation  (37)  the  frequencies  for  a  non- spinning  missile  differ  only  in 
sign,  (0^  »  -0g  «  0).  The  quasi -linear  solution  of  Equation  (52)  is 


a  special  fora  of  an  epicycle, 

8  " 

.  8ffi  cos  (0"p  ♦  03), 


101  102  _iB 
K  (e  1  +  e  2)e  19 

.1 


(53) 


i 

i 

i 

) 
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- ■■^...1^.«ky  .. 


■WtV'i  ■  ; 


for  which  Equation  (38)  reduces  to 


f 


“  Mq  “  ¥  60  «2 


(54) 


where  8 -  2K.  (maximum  value  of  8) 


K  m  ha  -  6  -  ®  - 


20* 


The  exact  solution  of  Equation  (32)  is  an  elliptic  function  with 
period  determined  hy  complete  elliptic  integral  of  the  first  type.  Thus 
by  use  of  a  little  algebra  and  a  table  uf  complete  elliptic  integrals 
it  is  possible  to  compare  the  periods  predicted  by  Equation  (54)  with 
the  period  of  the  exact  solution  of  Equation  (52).  In  doing  this  two  cases 
must  be  considered: 

1.  A  moment  vhich  grows  faster  than  a  linear  moment  (J^<£0) 

2.  A  moment  vhich  grows  slower  than  a  linear  moment  (Hf>0)  and, 

l“2  8  I 

therefore,  actually  changes  sign  for  I  -jjr-— I  ■  1. 


Although  a  quasi-line ar  approximation  may  be  reasonably  good  for 


reasonably  large  values  of 


cannot  be  good  for  values  of 


when  1^  is  negative,  it  certainly 
near  unity  when  is  positive. 


With  this  in  ilnd  we  may  now  state  the  surprln  -g  r.  u?  t*  of  this 
comparison  with  tivi  exact  theory.  For  negative  Kg's,  Equation  (82) 
predicts  the  period  with  less  than  1$  error  whenever  the  non-lineup  ^2 

K~ 

For  poeltlve  Jt,'»  the  error  will  be  less  than  a  percent  when  the  non¬ 
linear  moatent  contribution  is  lass  than  one-half  the  linear  lament 

for  larger  anglta.  As  a  result  of  these  facts  it  is  reasonable  to  make 
use  of  Zquatlomi  (3$  -  38)  with  considerable  optimism. 


■m 


i-:’ 


t 


?*<  .... 


Turning  now  to  the  question  cf  more  general  non-linearities  an 

examination  of  the  algebra  used  in  Equations  (55  -  58}  indicates  that 

2 

any  polynomial  function  of  6  could  be  used.  In  order  to  obtain 

2 

the  effective  values  of  higher  powers  of  6  we  look  for  those  terms 

i0l 

2n 

.  in  8  which  upon  multiplication  by  X  yield  terms  in  K^q  e  or 

^2 

K_-  e  .  Clearly  the  only  terms  which  have  this  property  are 

*u  i(s)  -  i  )  t(&  -Ci 

constants,  e  '^2  'l  ,  and  e  ,f±  rZ ' .  In  the  binomial  expansion 
2d 

of  8"  the  only  term  m  which  these  exponentials  appear  are  various 

r  1(0  -  02)  i(02  -  0rt)l 

multiples  of  the  cosine,  i.e.  je  +  e  I  . 

Although  the  constant  terms  are  unaltered  in  all  three  forms  of  effective 
yaws,  [62n]el,  P]e2  and  (p^,  the  cosine  ter®  makes  different  contri¬ 
butions  to  these  forms.  For  B2n  .  its  contribution  is  CK--2,  for  ftS2”!  - 

■  2  r2;n  k,  2  0.  -  **  0-  L  >2 

it  is  CK  2,  and  for  |62nJ  it  is  C  (— - i - - -). 

10  * 


As  an  example  of  this  algorithm  we  calculate  the  three  effective 


forms  of  8. 

Fro®  Equation  (25), 


i(0x  -  02)  i(0g  -  0J 


6*«  I^a  +  b(e  ^  +e  ■  *  ^  )J  (55) 

where  a  -  K^2  +  K^2 
b  "  *1*0  ‘ 

The  exponential  form  of  the  cosine  allows  an  easy  selection  of  the 
constant  and  cosine  terms  of  the  expansion.  These  particulT  terms  of 

I- 

the  expansion  will  be  identified  by  brackets  on  8* 

W  -  a2  +  2b2  +  2a  !  1  (56) 


M  -  a2  +  2b2  +  2a 


In  order  to  obtain  [B4]^,  [B4Je2,  or  \b\,  X,/.  K-Q2,  rx 
*10  “  *20  ^2 

- l - 1 -  respectively  should  be  inserted  in  the  empty  brackets  of 

Equation  (56). 


.  •  .  .  -t: 


^  V  ^  « 


;:>V?  ;  \  \  - 


,\i  /■'  • 

\  ’ 


[«  J.a  =  a2  +  2b2  +  2a  K^2]  (58) 

[s4]e  -  a2  +  2b2  +  2a  -1  '  ?f  ^  (59) 

L  h  -  4 

2q 

In  Table  I  values  of  8  are  given  for  all  valuea  of  n  between 
one  and  eight, 

4.  GEOMETRIC  NON-LINEARITIES  OF  LARGE  YAWING  MOTION 


The  basic  quantities  measured  In  free  flight  are  the  co-ordinates 
of  missile's  center  of  mass  and  the  direction  cosines  of  Its  axis  of 
symmetry.  Since  the  usual  formulas^  relating  the  direction  cosines 
and  the  motion  of  the  c.m.  to  the  conplex  yaw  X  were  derived  for  small 
yaws,  the  exact  relations  have  to  be  derived  for  this  report.  In  this 
derivation  we  will  find  it  convenient  to  keep  the  2-axis  in  the  horizontal 
plane  (fixed-plane  co-ordinates).  Unfortunately  these  co-ordinates  are 
not  the  some  as  the  non-rolling  co-ordinates  (v  «  0)  and  it  is,  therefore, 
necessary  to  calculate  $  for  these  co-ordinates.  Surprisingly  enough  in 
Reference  8  it  is  shown  that  not  only  is  v  finite  but  it  has  a  non-zero 
average  value.  In  our  development  the  cumbersome  Euler ian  angles  of  that 
report  will  not  he  used  and  the  desired  results  will  be  obtained  in  a 
somewhat  simpler  fashion. 

Co-ordinates  in  the  free  flight  range  system  will  be  Identified  by 
(x^,  Xg,  Xj)  and  in  the  fixed-plane  system  by  (y^,  y2,  y^).  The  range 
system  has  a  1-axis  pointed  downrange  along  the  intersection  of 
horizontal  plane  and  vertical  plane  containing  the  gun.  The  2-axis  lies 
in  the  horizontal  plane  pointing  to  the  left  and  the  5-axis  up.  The 
fixed-plane  co-ordinates  have  the  1-axls  along  the  miasile*s  axis,  the 
2- axis  in  the  horizontal  plane  pointing  to  the  right,  and  the  5- axis 
down.  Finally  the  non-spinning  co-ordinates,  which  ere  our  fundamental 
co-ordinates  in  the  theory,  have  the  1-axis  along  the  missile's  axis, 
the  2-arls  initially  pointing  to  the  right  in  the  horl.oncal  plane 
but  moving  so  that  $  la  zero  and  the  5-axis  fixed  by  tbe  rlgnx-L~cd 
rule. 


In  order  to  calculate  the  complex  yaw  it  is  necessary  to  know  the 
components  in  the  range  system  of  the  unit  vectors,  ?v  *3,  along  the 
fixed-plane  axes.  Since  the  fixed  plane  1-axis  lies  along  the  missile's 
axis  of  symmetry,  components  of  the  unit  vector  2  are  (n^,  n2,  n?)  where 
the  n,  are  direction  cosines  of  the  missile's  axis  with  respect  to  the 
range'system.  The  restriction  of  the  fixed  plane  2-axis  to  the  horizontal 
plane  is  equivalent  to  the  requirement  that  the  third  component  of  ^ 
be  zero.  This  together  with  the  requirement  that  it  be  perpendicular 
to  the  1-axis  and  pointing  to  the  right  completely  determine  the 
coatponents  of  the  unit  vector  along  the  2-axis. 


fa l  +  “2  V“l2  +  42 


According  to  the  right  hand  rule,  the  third  unit  vector  lr  equal  to  the 
cross  product  of  the  first  two. 

-a  v  #  °2  _ .  n?; _ 


*5  *  x  %  “  V  a3^  *  ( 


/  2,2 
K  +B 2 


-n/nl2  +  “22 


ifa.  +n2 


With  this  information  the  matrix  aquation  for  changing  from  the  range 
coordinates  to  fixed-plane  co-ordinates  can  now  be  written: 


rg  2 

“l  +  n2 


rr/r 

nl  +d2 


2  2 
+n2 


/2 2 
Pi  +a2 


2  2 
1  +n2 


;  iv  am.  A/g  idtrutnuiki iAiak 


M&iu. 


e£k**X±- 


If  the  vector  (±^,  x^,  x^)  is  inserted  on  the  right  side  of  Equation  (63), 
the  left,  side  would  yield  the  components  of  the  velocity  vector  in  the 
fixed-plane  co-ordinates.  These  co-ordinates  will  be  denoted  hy  the 
Symbols  (u^,  Og,  6^)  and  the  symbols  without  toe  circumflex  will  be 
reserved  for  the  non-rolling  co-ordinate  system. 


(64) 

(65) 


a  n2*l  -  nl*g 

**  f~~2  ",  T 

V11!  +  °2 

a  Wl  +  *2*3*2  *  (nl2  +  *2S)*5 

u,  ■ . . .  . . 

V11!2  +  *£ 


Dividing  by  the  magnitude  of  the  velocity  vector,  u  «  x^l  1+fa— 
/dxc\  2ll/2 

(  J,  and  simplifying,  we  can  obtain  the  following  expressions 

for  the  components  of  the  complex  yaw  in  the  fixed-plane  system: 

a*a- 

*  HU5i _  «*> 

(A2  ♦  °22>  Vi  *  <^r  *  <4?r 
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*2  “T 


_  _  <>t 
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For  the  case  of  a  flat  trajectory,  the  derivatives  are  the  same 
order  of  wgnituds  as  n^  and  n^  or  smaller  and,  if  we  make  use  of  the 

fact  that  n^2  ■  1  -  n^2  -  n^2,  we  can  obtain  the  first  order  approximation 
of  the  yaw  components 

A 

*2  -n2  " 


<*2 

*5 


A 

*3 
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Turning  now  to  the  problem  of  computing  the  axial  angular  velocity 
of  our  fixed-plane  co-ordinates,  we  first  must  state  the  definition  of 
toe  angular  velocity  vector  in  terms  of  unit  vectors  along  these  co¬ 
ordinate  axes  e^,  e£,,  ey 


where  (  )'  denotes  . 


Substituting  in  Equation  (70)  the  co-ordinates  of  the  vectors 
and  from  Equations  (61-62)  we  have 


(75) 


A  A  , 

V  *  «  ( 


7^; 


If  the  angle  between  the  missile's  axis  and  a  vertical  axis  is 
denoted  by  then  nz  »  cos  Y  and  Equation  (74)  reduces  to 

-A 

v  -  -  M.J  cot  ♦  (75) 

Since  Jig  is  the  angular  velocity  of  the  missile's  axis  about  the  2-axis 
which  is  fixed  in  the  horizontal  plane,  its  integral  should  be  related 
to  +•  Integrating  Mg  as  given  by  Equation  (73)#  we  see  that  this  is  true. 


Mg  dp  »  arccos  n^(p)  -  arecos  nj(o) 
o 

***  *  *  dp  *  * o 


Although  Equations  (75)  and  (77)  are  very  useful  for  the  numerical 
integration  of  the  complete  equations  of  motion,  we  will  find  Equation 
(72)  to  be  much  more  oeumiient  for  the  purposes  of  this  report. 

If  the  flat  trajectory  s-  -roxlmation  ■  pd  is  used,  a  first  order 
approximation  of  v  may  now  be  obtained  from  Equations  (67),  (68),  and 
(72). 

*  “  n 3n^  -  d3  +  t|)  (tg  +  -^) 


where  T 


A  good  approximation  to  the  yaw  components  is  an  epicycle. 
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where  K10  are  constants  and 


0^  are  linear  functions*  of  p. 

.  • .  i£  -  (%0  •!»  \  *  «m  •*»  \><-  *10  K  ““  *1  '  *»  ^  "" 

■  -  I  Ro*  K  *  &]  *  s  ["lO2  ^  “*  ^  *  **>  ^“4 

-  <K*^lK10lt20,1“^*ln^  <Sl) 

,OT  *,»«<»  (81)  It  <=«  f  •««»  that  v  h*a  *  «■-“«  «“•*  '*1"* 

<*  *  8  [*»  K  *  *»  5l]  '  (?ort»»«ly.  for  —H  J—.  <*•  f 

be  neglected.) 

Taming  to  the  wall' trajectory  term  T  ve  consider  only  the  effect  of 
gravity  and  lift  force  and  obtain  the  following  relations  for  constant 

velocity  uQ  from  Reference  5- 

*‘”\l 

j-vV-{^“‘,'X  co,\! 

where  are  seal  1  constant**)  and 

g  is  the  acceleration  due  to  gravity 

A  little  algebraic  manipulation  results  In  lie  relation  that 

*  The  circumflex  on  the  flR'a  Is  to  Indicate  that  the.,*  angles  are 
aaaaured  In  the  fixed  plane  system  and  not  in  the  non-rolling 
system  of  Section  2. 

**Thls  la  due  to  tbs  a« sumption  of  a  flat  trajectory. 
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+  fluctuating  terms  (84) 
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Since  for  most  missiles  10  ,  T  is  effectively  the  same  order 


0, 


as  and  we  need  only  consider  the  fourth  order  terms  which  are 
independent  of  the  trajectory  to  obtain  a  complete  fourth  order 
approximation  to  the  non- fluctuating  or  "D.  C.  component”  of  v.  This 
expression  is  calculated  in  Appendix  B.#  If  we  add  T  to  this  and  define 
f-yj  to  be  the  "D.  C.  component"  of  v,  we  see  that 
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With  these  results  we  can  now  calculate  either  the  exact  position 
or  the  average  position  of  2  and  3- axes  of  the  non-rolling  co-ordinates. 

Since  our  fixed  plane  co-ordinates  are  turning  at  an  angular  velocity  of 
^  with  respect  to  the  non-rolling  co-ordinates,  we  have  that 
P 

if  v  dp 

Lie  °  (86) 

As  it  ahall  be  shown  later,  v,  in  Equation  (86)  may  be  replaced  by  |vj  for 
the  data  reduction  of  large  yaw  firings. 

The  effect  of  large  yawing  motions  on  the  results  of  Section  £  manifests  it¬ 
self  in  two  ways.  First  the  distinction  between  the  non-spinning  co-ordinates 
and  the  fixed-plane  co-ordinates  becomes  important .  Accord!  'g  to  Equations 
(85-86)  these  co-ordinate  Bystems  rotate  with  respect  to  each  other  with  a  non-zero. 


*  In  Appendix  B  the  freqwncy  for  pure  precesslor  is  compared  with  the 
approximate  relation  baaed  on  this  expression  ior  rOl  .  As  can  be  seen 
in  Fig.  15,  the  agreement  is  good  for  8  *  sin  45°,  ■* 
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i  -  (|)  (52} 

Inserting  Equation  (25)  in  Equation  (90)  we  have 


1  "  ‘  5  *10  K20 


r  .  1(0, .  0 )  K0  -  0  >i 

(0[  -  02)  e  -  +  (0'  -  0X)  e  2 


Substituting  Equation  (91)  in  Equation  (89),  and  operating  on  Equation 
(89)  in  the  same  way  as  we  treated  Equation  (30),  the  fpllowing  two^ 
equations  may  be  obtained  from  the  coefficients  of  K^e  1  and  2 

respectively  0* 

<2-  V0’  +  M0  +  ^  4  t  0^)(-|)  -  0  (92) 


022-  v^+M0+4  Be2  +  K102(4  '  "  0  -  (95) 

Eliminating  first  Mq  and  then  "v  between  Equations  (92)  and  (93), 

4  *  *2  *5  fioX  *  «hoVT]  ■  7  *  fig]  |°->  <9‘> 
<  -4*1  -"o*^  ».*  <95) 

Thus  the  geometric  non-linearity  of  Equation  (89)  introduce  correction 

terms  to  the  left  sides  of  Equations  (57  -  38)  of  the  previous  section 
*  1 

and  replaces  It,  by  -  g  Mq  in  those  equations.  If  the  frequencies 

in  Equations  (94  -  95)  are  replaced  by  their  fixed  plane  values  from 
Equation  (87  -  88),  and  fourth  powers  of  K^'s  are  neglected.  Equations 
(94  -  95)  may  be  written  in  the  following  useful  form. 

'  >^(-^'>°g  ->  <(«) 
L  J  0i  -  h 

K  •  K  |i  -  5  (^O2  +  j  *  Mo  1  6e2  (97) 
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An  examination  o f  Equation  (97)  shows  that  the  effect  of  a  large 
percentage  change  In  the  ftequmlM  Implied  V  Equations  (87  -  88)  Is 
almost  completely  cancelled  vhen  the  terms  In  the  cosine  of  are 
properly  handled.  A.  a  Chech  on  the  derivation  of  Equations  (96  -  97)  they 
are  derived  in  Appendix  C  directly  from  the  yaw  equation  expressed 
In  the  fixed- plane  system.  This  derivation  show,  that  the  assumption 
implied  in  Equations  (67-83)  Is  equivalent  to  the  basic  assumption 
underlying  this  whole  report. 

5-  EXPERIMF-WTAu.  results 

In  Section  3,  a  wry  simple  example  or  non-linear  yawing  motion  with 
two  degrees  of  freedom  Is  considered  in  great  detail.  The  very  same 
quasi-line ar  assumption  may  be  applied  to  a  much  more  complicated  system. 
In  appendices  D  and  E  the  yawing  and  swerving  motion  of  a  missile  for 
which  til  aerodynamic  coefficients  are  quadratic  functions  of  t  la 
considered.  Equations  relating  the  parameters  of  the  epicycle 
approximation  (frequencies  and  damping  exponents)  to  the  coefficients 
of  tbs  parent  non -linear  equation  and  the  amplitudes  of  both  modes  are 
derived.  Since  these  relatione  are  essentially  approximations,  their 
value  must  be  determined  by  either  numerical  or  experimental  checks. 

A  good  check  on  these  approximations  can  be  obtained  by  means  of 
the  actual  free  flight  motion  of  a  mloella  acted  on  by  non-linear  forces 
and  momenta,  for  this  reason  various  results  of  firing  taste  on  «'• 
.park  range.9  were  earned  for  possible  verification  of  the  quael-llnoar 
hypothesis.  The  result,  of  this  Investigation  eerw  quite  encouraging. 

The  affect  of  a  cubic  static  moment  on  three  different  programs  was 
first  considered.  These  program#  were  a  finned  missile  program  (81* 
gjTtsn10),  a  body  of  revolution  program  (Army-Eovy  flpinner  Rooxet 
program6),  and  a  large  yaw  body  of  revolution  program  flrdd  by  E. 

Koecker. 

in  the  finned  missile  program  a  free  flight  rente  value  of  was 
computed  from  the  product  of  the  frequencies.  According  to  Equation 
(36a)  and  the  definitions  of 

(56> 
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Thus  the  range  value  of  K^,  as  obtained  from  an  epicycle  fit  of 
Individual  rounds,  should  be  a  linear  function  of  the  corresponding 
effective  squared  yaws.  In  Figure  2  the  range  values  of  are 
plotted  against  &  ^  and  different  lines  for  each  configuration  tested 
are  dravn  through  corresponding  points.*  The  linearity  exhibited  by 
these  experimental  points  constitutes  the  first  check  of  the  theory. 

It  should  be  emphasized  that  In  all  cases,  epicycles  with  damping 

were  fitted  to  the  data  and  although  the  exponentially  varying  K. *8 

2  1 

actually  should  have  appeared  In  the  calculation  of  &e  they  were 
approximated  by  their  midrange  values  K^q'3 .  The  approximation  is 
Implied  by  our  basic  quasi-llnear  assumption. 

The  Arny-IIavy  Spinner  Rocket  program  consisted  of  three  model 
lengths  each  with  three  different  center  of  mass  positions  and  fired 
at  three  Mach  numbers.  A  maximum  of  twenty- seven  values  of  was 

therefore  possible.  Since  for  each  case  at  least  four  points  with  a 
2 

reasonable  qptsd.  In  are  necessary  for  a  good  determination  of 

it  was  actually  possible  to  obtain  only  sixteen  values  from  a  careful 
analysis  of  the  126  record  rounds  of  "he  program. 

These  values  vith  their  standard  errors  are  listed  In  Table  II. 

A  sample  plot  of  the  experimental  points  for  the  9-callber  model  at 
Mach  number  of  1.8  la  given  in  Figure  3* 

For  some  of  the  Mach  numbers  and  configurations,  wind  tunnel  data 
taken  by  R.  Krleger  w ere  available.  Hie  data  for  the  overturning 
moment  were  fitted  by  cublcs  and  the  cubic  coefficients  (K,,  2)  axe 
listed  in  Table  II.  At  all  seven  points  of  comparison  the  agreement 
Is  good. 


Since  these  finned  missiles  were  not  spinning,  0*  •  -  and  the 
effective  squared  yaw  assumes  the  concise  f on  ^  (K^q2  +  K^Q2. ) 

In  Reference  10,  the  Kj,  ’s  were  erroneously  plotted  against 
"range 

“5  2  2 

6  "  K10  +  Kqq  and,  hence,  the  slopes  obtained  In  thai,  report  are 


actually 


2* 
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TABLE  II 


Mach 

Rusher 


1.3 

1.8 


1.3 

1.8 

2.5 


1.3 

1.8 

2.5 


FOR  THE  ARMY-HAVY  SPIHNER  ROCKET 


Forward 
Center  of  Mas 8 

Range  Wind 
Tunnel 


Middle 

Center  of  Mass 

Range  Wind 
Tunnel 


Fl,-e  Caliber  Long  Models 


-  2 
-14 


3  +  1  0 

-7 


Seven  Caliber  Long  Models 


.21+1 
-12  +  1 
-28  +  1 


-11 

-26 


-  5  +  1 


wine  Caliber  Long  Models 


-  35  *  1 
-30  +  4 


—  16  +  1 

-  30.+  2 


Rear 

Center  of  Mane 

Range  Wind 
Tunnel 


1+2  3 

3  +  1  0 


7+2  5 

0  +  3  5 


.1  +  .01  * 

-  5-7  ♦  -5 

-  1  +  .6 


*  Tt*  c-ttr  «  ««•  BOK  mUU  U  .2  «1.  r.«  irf  «• 

XOCfttlOP* 


tsB&ssae^*s#t 


Since  models  with  throe  different.  center  of  mass  locations  were 
fired,  a  second  check  is  possible.  If  the  normal  force  is  expanded 
as  a  cubic  function  of  6,  the  usual  center  of  mass  relations  provide 


V°  ■  \  * 4  *»„  <99) 

*,.,2^  ’  *S^2  *  11  %62  (100) 

where  Kjj  (4)  are  the  moment  coefficients  for  a  center  of  masr  located  q 
calibers  recur  of  the  center  of  mass  for  the  K.,  *s. 


This  means  that  the  K.,  and  K,,  8  for  different  center  of  mass 


\“4V 


locations  are  linear  functions  of  location.  It  was  possible  to  measure 
three  values  of  K„  2,  for  on-ly  one  configuration  and  one  Mach  number. 


In  Figure  4  both  and  5C,  _  for  this  case,  are  plotted  against  center 
”0  V 

of  mass  location.  The  fact  that  the  cubic  coefficients  as  well  as  the 
linear  coefficients  fall  on  a  straight  line  is  another  point  in  favour 
of  the  theory.  The  slopes  of  these  lines  are  Kjj  2  and  respectively. 

For  six  other  cases  it  was  possible  to  compute  2  from  two  values 
of  Kj,  2  and  all  seven  values  together  with  the  corresponding  linear 


values  (jc.  )  are  tabulated  in  Table  III.  Once  again  Krieger’s  wind 
0 

tunnel  data  were  analysed  and  it  was  found  that  at  all  points  of 
comparison  the  wind  tunnel  results  were  In  good  agreement  with  center 
of  mass  values. 

In  Appendix  D,  it  is  shown  that  the  lift  force  and  Magnus  force 
coefficients  may  be  directly  measured  from  the  swerving  motion  when  it  is 
large  enough.  According  to  this  appendix,  most  of  the  swerve  is  associ¬ 
ated  with  the  lower  frequency  and,  hence,  the  rat^e  values  of  and  K^, 

2 

for  each  round  should  be  plotted  against  the  corresponding  value  of  &e2  • 
A  sample  plot  of  this  process  is  given  in  Figure  5- 


42 


it  is  necessary  to  expand 


In  order  to  calculate  g  from  K^g, 

Equation  (17 )  in  powers  of  62  and  compare  coefficients 

vv\ 


v-v*  w 


Using  these  relations  it  was  possible  to  calculate  eight  values  of 
both  and  g  from  the  swerving  motion  and  these  are  given  in  Table 

III.  With  the  exception  of  two  values  of  g  for  the  7-caliber  models 

the  agreement  is  remarkable.  The  reason  for  these  two  discrepancies 
is  at  present  unknown. 

Although  the  lateral  displacement  due  to  Magnus  force  is  about  one 
tenth  that  due  to  lift,  it  was  possible  to  make  two  measurements  of 
Kp  g.  The  experimental  points  for  these  two  cases  are  given  in  Figure  6. 
ThS  excellent  Internal  agreement  of  these  quite  delicate  measurements 
Is  extremely  gratifying. 

Finally  in  Roecker's  large  yaw  program  It  wan  possible  to  check  the 
treatment  for  the  geometric  non-linearities.  In  Figure  7  KL,  ,  as 


calculated  from  Equation  (97),  1*  plotted  against  6  .  Here  we  see 

that  the  data  la  essentially  bilinear.  Each  line  corresponds  co  a 

cubic  segment  in  the  moment  plane.  If  the  parameters  of  each  cubic 

are  calculated  from  the  slope  and  intercept  of  its  corresponding  line 
2 

in  the  plane,  they  can  be  pieced  together  to  for  t  a  smooth 

moment  plot.  (Figure  8).  (An  examination  of  the  spark  shadowgraphs 
revealed  the  fact  that  flow  separation  occurs  at  about  21°  and  this 

explains  the  sudden  change  in  the  moment  curve  it  this  point.)  In 

o 

Figure  9>  is  plotted  versus  Beg  and  the  corresponding  lift  force 
plot  is  given  in  Figure  10. 


Since  wind  tunnel  measurements  for  this  conflgur*!!?”  had  been 
made  by  V.  Buford,  a  comparison  was  possible.  The  data  was  divided 


Range  valuta  of  Kg  -  and  for  tba  Army- Nary  Spinner  Rocket 

program  were  fitted  by  least  squares  to  Equations  ( 10 5- 1(A)  and  values 

* 

of  Kjq,  KgQ  -  K^,  Hj,  2  11114  Kg^2  ver®  C0B3?u"te<1,  111  all  cases  the 

coefficients  of  Kg^g  were  small  and  this  quantity  was  poorly  determined. 
K^2  v&b  therefore  oelttcd  from  Equations  (103-104)  and  they  were  fitted 
separately  to  rac^t  valuta  of  Kg  -  and 

Since  this  ■eatureawt  depends  on  the  damping  exponents  and  hence 
la  qalte  delicate.  It  vaa  possible  to  obtain  only  eight  values  of  Kp  2 

from  K_  and  four  from  (K_  -  £,.)  .  These  are  listed  in  Trbi*  V 

range  range 


into  angles  of  yav  lens  than  21°  and  greater  than  21°  and  pairs  of  cubics 
fitted.  The  resulting  coefficients  are  compared  in  Table  IV.  The 
agreement  is  excellent. 

As  a  result  of  this  rather  spectacular  success  with  the  non-linear 
static  moment,  the  more  difficult  problem  of  Magnus  and  damping  moments 
was  then  considered.  In  Appendix  E  it  is  shown  that  a  quadratic  dependence 
of  Kj,  and  Kg  -  on  8  so  affects  the  damping  of  the  epicycle  that  the 
usual  linear  formulas  for  Kg  -  and  are  acti:ally  the  following 
combinations  of  coefficients: 


(*«-  *»> 


Sim  *  \2  <  f 


-  v  $  (|f |>  <ki°2  -  ^ 

,  .  a  . y.  \oK2  -  *x,%2. 

*TO  \2  6e  *  Vl}  (  jj.  2  .  ^  } 


W  2. 

*  *^2  "  Sttg2  +  ?  KMAO 


TABLE  IV 


Large  Yaw  Results  for  M  =  2.3 


V 

5k 

V 

\ 

V 

Wind  Tunnel 

1.85 

.625 

-  1.1 

1.01 

L.O 

■before  separation 
f  x  „  0 

Range 

1.82 

.829 

-  1.0 

0.99 

3-6. 

|  8  ?  sin  21 

Wind  Tunnel 

.129* 

0.0 

1A5* 

0.0  ' 

after  separation 

Range 

.719* 

0.1 

1A3* 

o.k  J 

8  «  sin  21° 

*  Eveduated  at  separation 


TABLE  V 


Mach 

Ho. 

Ky  2 

f5 

V 

<V  V 

JV 

C.M. 

Swerve 

FCM 

MCM 

RCM 

Seven 

Caliber  Long  Models 

1.3 

lit 

1?  +  1 

lt.lt 

1C  -  16  +  5 

range 

•  5  +  1 

OC  -  1C.)  -20+14 

“  ™  range 

— 

— 

1.8 

22 

— 

4.8 

Km  -34+16 

range 

— 

0 

25  * l' 

— 

— 

2.5 

19 

— 

5-1 

K.  -36*5 

range 

— 

-6 

(1C  -  1C,)  -28  +  4 

“  m  range 

— 

•** 

Fine 

Caliber  long  Models 

1.8 

12 

6.3 

1C  -29+2 

range 

— 

-5.2 

0C  -  1C*)  -26  +  3 

«»■» 

range 


r  p  for  Mach  number  1.3  vns  given  In  Figure  6  as  8  +  1. 
'a2 
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SOI.  Wind  Tunnel  meaaureaento  were  made  by  Uicbuk  and  Spt.  rfci'  f1 or 7 
caliber  model  at  M  between  1.6  and  2.5*  According '-to  tht.r  vaults 


and  the  agreement  between  pairs  of  K^,  2  is  within  the  standard  errors.*- ■ 

6 

In  Figures  11  -  12,  a  pair  of  plots  of  experimental  points  for  K_ 

1  range 

and  (K_  -  K,.)  are  given, 

range 

In  all  four  cases  values  cf  K_,  „  for  two  center  of  mass  locatio.  s 

A6* 

were  obtained  and  K^,  £  s  and  CP_  2'  B  were  computed  from  these.  At 
8  r8 

the  Naval  Ordnance  Laboratory,  Luchuk  and  Sparks11  have  made  wind  tunnel 
measurements  of  these  quantities  and  their  results,  as  given  at  the  bottom 
of  Table  V,  are  in  reasonable  agreement  with  flight  tests.  Finally  l-; 
should  be  noted  in  the  Table  that  at  one  point  of  comparison  of  cental  of 
mass  and  swerve  values  of  K^,  2,  (  the  seven  caliber  models  at  M  ■  1.3.  the 

agreement  is  good. 

As  a  last  example  of  this  quasi-linear  technique  we  will  consider  i 
body  of  revolution  which  displayed  an  extremely  non-linear  Magnus  warn'-. 
The  range  values  of  for  this  model  were  obtained  by  E.  Roecker  and  are 
plotted  versus  8  2  In  Figure  13.  Although  these  data  are  fitted  by  two 
lines,  only  the  first  line  is  well  determined.  In  fact  if  the  fourth  po’ at 
from  the  right  were  neglected,  tne  large  yaw  values  would  be  reasonably 
well  represented  by  the  dashed  horizontal  line. 

In  Figure  14  the  corresponding  cubic  segments  are  plotted  and  compared 
vtth  a  BRL  wind  tunnel  curve  obtained  by  A.  Platou.  The  good  qualitative 
agreement  for  this  strongly  non-linear  moment  Is  remarkable. 

Since  the  quasi-linear  relatione  have  been  so  well  ve  'tf led  by 
experiment,  they  should  be  quite  valuable  In  the  Important  problem  of  the 
prediction  of  missile  motion  from  a  knowledge  of  the  force  and  moment 
curves.  This  application  of  the  quasi-linear  technique  is  described  In 
more  detail  In  the  next  section. 


Thus  an  observed  dependence  of  (Kg  —  K^)  or  ■’wjciltude  of  yaw  le 

range 

not  necessarily  due  to  non-linear  moments  but  may  be  due  to  a 

non-linear  Magnus  moment. 
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6.  PREDICTION  OF  NGN- LINEAR  YAWING  MOTION 

Although  the  application  of  the  theory  of  this  report  to  range  data 
analysis  enhances  markedly  the  value  of  spark  ranges,  an  even  more  important 
application  is  the  prediction  of  the  yawing  motion  of  a  missile  acted  on 
by  non-linear  forces  and  moments.  The  experimental  results  of  this  report 
indicate  that  this  inverse  problem  should  be  successfully  handled  by  the 
same  theory.  In  this  section  we  will  outline  the  procedure  for  this 
prediction  problem. 

First  the  sizes  of  the  two  arms  must  be  obtained  from  the  initial 


conditions 


**•  N>“*l 


,  ,  10.  ,  10p 

X0  -  (-  *  10^  e  X  +  (-  +  i02)K2  e  *  (106) 

In  order  to  take  care  of  the  effect  of  damping  the  trajectory  should 
tie  divided  into  intervals  over  which  neither  amplitude  changes  by  more 
than  5<#.  The  length  of  the  first  interval  may  he  estimated  from  the 
linear  draping.  If  the  calculated  quasl-linear  damping  is  much  larger, 
the  interval  may  then  be  shortened  and  the  process  iterated.  The  lengths 
of  the  other  Intervals  should  be  determined  in  a  similar  manner.  Values 
of  0  ,  frequencies  in  non- rolling  co-ordinates,  can  now  be  computed  from 
Equations  (92  -  93)  which  are  simple  quadratic  equations.*  may  be 
computed,  from  Equations  {87  -  88).)  Then  value  si  of  may  be  computed 
from  Equations  (e4  -  E5)  which  are  even  simpler  linear  equations.* 


*  For  ease  of  calculation  the  linearized  values  of  P,s  and  a.s  may  be  first 
placed  in  the  coupling  terms  arising  from  the  ncm-Ilnear  geometry  and  an 
Iteration  performed  If  needed. 


*  I  « 


If  a  more  complex  polynomial  dependence  of  and  K^,  on  &  is 
required,  this  may  he  handled  in  the  following  way. 


%  -  y v*  ;io6) 

k  =  0 

Then  Equations  (92  -  93)  and  <S4  -  E5)  take  on  the  following  slightly 

* 

more  complex  form: 

n  S 

y.  4  <>*'*  *  *»  <4  ■  4><r>  *  “*(H  •  “i'  * 0  (109> 

k  -  0 

.o  \  n  #  *  ,«2kx  .  _  ri5  wlh  +  «,(#  -  Co)  a  0  Clio) 


012-  V02 


(62k)e2  ♦  Kx02(^2  -  0X)(r)  +  °2(U  *  1 


-  »>  -  vi  y  4<»a,>.1  -  ^  d**  *  *»**'  * 


-  '  0^K2O2  +  ^l^^O2  +  (U1) 

.  - « 

Oj(2(lj  -  .)  •  »  y  4<6a».2  -  V  [<>'  ‘ 

*  Eqa.  (92-93)  are  Use  slightly  modified  by  the  present  of  damping. 
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and  those  terms  appearing  In  the  right  hand  summations  which  have  powerB 
of  8  greater  than  2n  are  neglected.* 

It  should  he  remembered  that  these  equations  are  derived  on  the 
assumption  of  a  flat  trajectory.  The  complications  introduced  by  gravity 
vlll  be  considered  in  a  later  report.  In  any  event  for  a  large  class  of 
problems  this  simple  procedure  should  allow  reasonably  good  prediction  of 
the  motion  of  mlaallea  acted  by  non-linear  forces. 

7-  SUMMARY 

1.  A  convenient  expansion  of  the  aerodynamic  force  system  for  large 
yaw  has  been  obtained  and  the  necessary  geometric  relations  for 
these  large  yaws  derived. 

2.  Relations  for  parameters  of  a  linear  equation  which  la  a {divalent 
to  the  actual'  non-linear  equation  have  been  derived. 

3*  These  relations  have  been  tested  by  actual  firing  tests  and 
excellent  Internal  consistency  has  been  observed.  Where  wind 
tunnel  measurements  were  available,  good  agreement  has  been 
obtained. 

*  ‘2 
The  effect  of  the  terms  in  JT  which  are  higher  order  than  8  on  the 

t  " 

A  term  In  Equation  (20)  has  been  neglected. 
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Jt  ?  ($2  '  ^1^K10  +  +  “2^0  ^  l'112^ 


2 


4 e2k  ■  4- 


*(JL2k  '  kl  JT.2k*  8 


4  *6*  lppl««lcc  of  thi.  ttcbaM*.  «  «*  ',f  th* 

,  mlBRlle  -ted  on  by  non-linear  forces  and 
ys.vd.ns  motion  of  a  missile  a-xea  on  j 


moments  has  been  described. 


Am.NPlX  *  :  Derivat  ion  of  Equations  of  Yawing  Motion 


In  this  appendix  the  exact  equations  of  yt.vi.ng  motion  vi  1 1  be  derived 
In  two  steps.  First  the  ,ir«s  equation  which  governs  the  magnitude  of 
the  velocity  vector  will  be  obtained  and  then  the  equations  for  the 
trr-nsveree  velocities  and  angular  velocities  components  will  be  derived. 

3y  means  of  these  equations  the  exact  equation  for  the  complex  yaw, 

X,  may  be  formulated. 

The  fundamental  force  equation  for  the  velocity  vector  may  be 
written  In  the  forms 


(«) 


where  ?  is  the  aerodynamic  force  vector, 

l?  is  the  acceleration  due  to  gravity  vector. 
If  th*  s  equation  la  dot  multiplied  by  u, 

—a  dlf  du  4  *4  4 

BW-xr-*»3t”u-r+a7'« 


31 

du 


(A2) 


where  is  the  derivative  of  the  scalar,  u. 


Dividing  Equation  (A2 )  by  — j-,  changing  to  the  Independent  variable  p 

^  dt,  and  denoting  the  trajectory  component!  of  t  and  $  by  and  gj, 
we  can  obtain  the  final  form  of  the  drag  equation. 

«'.  V  .  *T4 
5 


fM) 


'  JD  *  Jg 


where  JD  -  Kjj 

j_  -  ^r4 

‘  T 
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Turning  nov  to  the  consideration  of  the  components  of  the 
velocity  and  angular  velocity  vectors  In  our  moving  co-ordinate  system, 
we  writs  the  equations  for  uhe  conservation  of  linear  and  angular 
momentum  In  the  fallowing  forms. 


L  (mt)  *t+ mi 


where  t,  •  K  angular  momentum  vector 


A  -  axial  moment  of  inertia 
B  -  transverse  moment  of  Inertia 
if  m  aerodynamic  moment. 

In  order  to  calculate  the  derivatives  In  Equations  (A4  -  A5)  for 
components  in  our  moving  co-ordinate  system,  a  useful  vector  Identity 
for  use  In  a  moving  co-ordinate  system  with  angular  velocity  if  will  he 
stated. 

It  @)  -  ^  ,  aFT>  +  ^  *  CA6) 

Applying  Equation  (A6)  to  Equations  (Ak  -  A5)  with  "5  ■  (J^,  av,.  «r  ), 

A  *5^,  M  M 

*  I  *2  *•}  “l  '  “l  “3  1  '  I  Tl  I  "  /  %  I  '*7) 


>  *2  ♦*,«!  - 

\  S  +  *i  *2  “  “l 

h  ^ 

+®3^1  '  *“3 

t  Bij  ♦  QjEtag  -  Aa^  J 


Since  the  axial  component  of  the  angular  momentum  equation  Is  uncoupled 
from  the  other  components,  we  consider  It  separately.  If  is  replaced 
by  Its  definition.  Equation  (7),  the  equation  for  the  axial  spin  may  be 
written  la  the  following  dimensionless  form. 


where  v 


axial  radius  of  gyration  In  calibers 


j  *  £i-  x 

A  ZB  A 

Substituting  In  Equation  (A9)  for  2  from  Equation  (A3),  we  can  write  the 

*  u' 

final  form  of  the  roll  equation. 

v  <■  (D  -  Jg)v  (AlO) 

_2 

where  D  «  -  k^ 

In  the  treatment  of  the  vector  components  in  the  plane  perpendicular 
to  the  missile' s  axis,  the  rotational  symmetry  of  the  force  system  may  be 
exploited  by  making  use  of  complex  numbers.  Multiplying  the  3- components 
by  1  and  adding  to  the  2-  components. 

Eg  +  IF.  (A11> 

(«2  ♦  lij)  -  ♦  i®;)  ♦  1^l(u2  +  iuj)  - - J -  +  ®2  +  ^3 


A®.  K,  ♦  1M, 

(«a  +  fc&j)  “  i(-g —  ^(©s  ♦  * - g - r 


In  tbs  analysis  of  this  report,  the  small  change  In  v  predicted  by 
Equation  (AlO)  will  be  neglected  and  all  derivations  will  assume 
constant  spin. 


(82+l85)4  t  , 

7  - - 5 

u 

^  ,  jj^  transverse  radius  of  gyration  in  calibers. 

fc  Keierenc.  3  it  i.  «-  *  -»  ««  '“'TltS 

*.  „  .OTr«ci.hl.  e«o=t  on  **  u»~is.d  ^  »«»  *» 

eete<ent  ,  .  Kr  thlo  nw  the  other  foroe  eoeftlclente 

Equation  (A15)  MU  »  negleoted  In  the  derivation  of 
Equation  (AW)  o»  no.  he  re-«rltte«  In  .  «oh  «re  oonveni 

Uu  .  \  *  («t  ♦  iv)i  •  ’  <U7) 

Differentiating  Equation  (A17)  end  Bolting  for  lAu* »  .  ^  , 

I,'.!-  1  X-.fc.rf-i'*] 

H  Equation  (U6)  1.  -It^U-  V  «  .  -  W  '  161  *” 

uaed  to  elialoato  “  **  *“ 

required  for  the  report. 

X"  +  |h  ♦  Jg  '  "1  “  1(7  "  2vJ  x 

#5 .  J)  .»(*»-«(*♦  V  i ’*  5’3  x  ("9) 

.  0  -  -  lv)7 

where  H  «■  -  JD  +  \ 


M  -  tej'2  JK  '  8jS  “  JD  +  V  "  **2  Jm 

T  *  ^ 
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APPENDIX  B:  lfth  Order  Approximation  of  v 

dx  dx 

If  the  trajectory  terms  are  neglected  (~ —  =  -  0),  Equations 

GX.  0^ 

(66  -  67)  for  the  two  components  of  yaw  in  the  fixed- plane  system 
reduce  to: 


*2  “  ~r 


2  2 
i  +  n2 


^[1  -  n32j) 


in3  L  f 1  n22  -  c32] 

-J-rVl 

L  1  •  n3  J 


Solving  Equation  (B2)  for  n^  by  an  application  of  Equation  (Bl) 

A  a  3  -  I/2 

n3  -  Xj  (1  -  Xg2)  (B3) 

A  1  A  A  g 

nj-Xj+^XjXg  +  .1.  (B4) 

Equation  (Bl)  is  now  solved  for  Hg  and  n^  replaced  by  its  value  from 

0*). 

A  1  A  A  9  ,  _ , 

®2  *  Ag  *  2  ^2  +  •••  (B5/ 

Row  from  Equation  '72)  we  know  that 
A  ~  ninln2^ 

“l  W  +  n2  ^ 

n3  fnl2n2  +  *  nxnj)n2^ 

*  - 7“S - 37"“ - 

\  (»!  +  n2  ) 

n3n2  n2a32n3 

•  — —  + - * - vBo; 

°1  n^2  +  ng2) 


mss 


|T'\  r-  ’7  :  r' ‘ \  S 

Is, 


'/■>  1/2  2X\  ,  2  » 

»  n^n^l  +  ^  (n,  +  nj  ;;  t  righ^  +  .  .  . 

Substituting  in  Equation  (B6)  from  Equations  (B4  -  B5), 

A  At  A  A  i  A  2  A  . 

v  «  Xj  +  Xjj  X2  Xj  +  .  .  .  ( 

As  in  the  text,  the  components  of  yaw  may  be  approximated  by  an 
epicycle. 

*2  «  *10  cos  +  XgQ  cos  $2  ( 

Xj  -  K^j  sin  ^  +  Kg0  sin  ( 


where  K^n  are  real  constants  and 

01Q  are  real  linear  functions  of  p. 

If  Equations  (S8  -  B9)  are  used,  &  small  algebraic  effort  puts 
Equation  (B7)  in  Its  final  form:  (BIO) 


v  *  -  (5) 


,  2  *•  „  ho2  *  ®a>#%  _  2{*  ,,  ^0  +  *10  v 
x0  ?!  U  + - 5 - >  +*20  ?2  (1  + - 5 - } 

mm  m 


+  periodic  terms  ♦  higher  order  terms 

In  order  to  get  an  estimate  of  the  range  of  accuracy  of  this 
approximation  we  will  consider  the  case  of  constant  precession  about  the 
trajectory.  Under  this  assumption  one  of  the  amplitudes  of  oscillation 
vanishes  and  for  a  linear  moment  Equations  (109  -  110 )  reduce  to 


*i'2  - 1  <V  *214  ho1  - 0 


*  2 
I  IS  6 


£4 


'Ll  *10*  *  “o 


‘•a  -  It,„a)l/2  ■  i  -  i  k  2  -  i  k,*  *  .  . 


From  Equations  (87  >68)  and  (BIO), 


K'  s  k'  ■  <I)K102  V  l1  *  t  Ku  > 


•  (B12) 


#l‘  AV  (i  'iKi02  -  JO’  (B13) 

i  ^  (1  *  |  KiQ2  +  |  K10U) 

Combining  Equations  (Bll)  ani  (B13), 

?i‘  -I  [l  ^l]  11  *  ®  *“  * *  K“)  <al,) 

la  taken  for  1  »  1  while  the  -  sign  is  appropriate  for  i  •  2. 

If  0^,'  is  the  rate  of  precession  about  the  trajectory  or  equivalently 
the  rate  measured  in  the  range  co-ordinate  system,  it  la  easy  to  see 
that  0^'  is  its  projection  on  the  Missile's  axis. 

•’*  ^ir'  -T^i')  <“-5) 

Although  v  is  actually  a  fluctuating  quantity  the  period  of  the  motion 
in  the  fixed-plane  co-ordinates  must  be  the  same  as  that  in  the  range 
co-ordinate  system.  Equation  (B15)  is,  therefore,  the  exact  form  of 
Equation  (B13).  As  can  be  easily  verified  the  expression  in  in 


Equation  (B13)  is  actually  an  expansion  of  I 


(1  -  ) 


If  the  first  and  second  order  forma*  of  (B13)  are  compared  with 
Equation  (Bl4),  a  measure  of  the  basic  accuracy  of  our  technique  as 
applied  to  the  geometric  non-linearities  is  possible.  Tn  Figure  15  the 
ratio  of  the  approximate  frequencies  to  the  exact  are  plotted  for  both 
fast  and  slow  frequencies  and  both  first  and  second  order  approximations. 
For  angles  up  to  45°  the  second  order  approximats  relations  show  good 
accuracy. 


■  tj  j*  -  f(l  -  5  K102)]  (1  +  2  Kl02)  while  the  second  order 


form  includes  terms  in  K 


m 


APPENDIX  C:  Effect  of  Large  X aw  In  Fixed- Plane  Co-ordinates 
If  the  overturning  moment  it>  the  only  aerodynamic  mnment  to  he 

2 

considered  and  this  moment's  coefficient  is  assumed  to  be  linear  in  5  , 
Equation  (A19)  assumes  the  form 

A."  -  +  i(v  -  2v)j  X  -  ^(Mq  +  &2  -  Q(v-v)  +  i(0  j  -  v  X  =  C 

For  the  fixed- plane  system  v  is  given  by  FiCuation  (8l^ .  Neglecting 
squares  of  Kj2  and  making  use  of  Equations  (90)  for  , 

X"  +  |^)(&2)'  -  i(v  -  2v^J  V  -  (Mq  +  V*  S2  -  vv  -  tv^  \  «  C.  (C2 

From  Equation  (78), 


*  (X  *  $  )&  -  $) 

-  -  •£  <xx  -TR'-V+lflf)  (C3 

If  w  consider  only  the  constant  terms  and  those  vlth  frequencies  (j^  -  $ 
and  -  ^)*  the  epicycle  solution  may  be  substituted  In  Equation  (C3) 
and  ve  obtain  the  following  result  n  a  *  .  (Ci 

A  l  T  22'  22'  .n.  +  ^2\  ”  ^1^71 

0  -  •  S  [*10  *1  +  *20  h  +  <-iT-2)K1K2  <e  +  «  )] 


2  -  ^  2)  i($x  -  ?2)  i(^2  -  J,) 

—  ■  "y-  *3*2  (e  *  • 


)  +  .  .  .  (C 


Substituting  the  epicycle  and  Equations  (C4  -  C5)  In  Equation  (C2)  ve  can 
obtain  a  pair  of  equations  similar  to  Equations  (92  -  93)  but  with 
additional  terms  coming  from  v  and  0  . 

K 2 ^  ».!2  <2;  -  K>  ■  ^  <=3;  -  vi 

ax ♦  s:2  *  i:. 


2 

A*  2  _flt  if  2  "lK10  .A* 

P2  -  v02  +  M0  +  Mg  fce2“  +  -M2-  <P0 


Pj.) - 5 -  (202  *  v^ 


K10 

r~ 


(C7) 


Equations  (C6  -  Cj)  may  now  be  simplified  by  use  of  the  linear  approxi¬ 
mation  v  .on  the  coefficients  of  the  quadratic  terms  K^q  . 

%(V  -  $[)  +  (|)K102  ^  -  $‘ )  -  Mq  +  l£  Bel2  (C8) 

%(y  -  -  (5)K202  K  (^1  "  “  *0  +  6e22  (C9) 


Eliminating  first  Mq  and  then  v  from  Equations  (C8)  and  (C9),  the  following  I 

w  i 

final  results  are  available. 


Since  these  equations  are  identical  with  Equations  (96-97)  the  v 
correction  for  the  co-ordinate  system  used  In  the  text  Is  based  on 
the  same  equivalent  linearization  assumption  as  the  other  results  of 
this  report. 

APPENDIX  D:  SWERVING  MOTION  FOR  NON-LINEAR  FORCES 

Since  a  proper  treatment  of  the  swerving  motion  should  yield  values 
0fKN,  2  which  may  be  coopered  with  the  results  of  the  center  of  mass 
methoS,  we  will  consider  In  the  appendix  the  swerving  motion  of  a 
missile  acted  on  by  non-linear  aerodynamic  forces.  If  the  aerodynamic 
force  in  the  range  co-ordinates  Is  denoted  by  (F^..,  Fgy,  F^r ),  it  can 
be  shown  frost  Equation  (63)  that 


61 


'1  -  -yl  -  n, 

'  \  R2  n3 


n5  0 


-  n.  F9  +  n-n^F,  / - 

F2r  ♦  “jr  '  <”2  *  *”3)ri  *  ,  ‘  M/*  '  *s‘  PJ 


Substituting  from  Kquations  (5-6)  and  neglecting  the  snail  1U  and 
terms  P  . 


T:' 2r  *  “Sr1  •  -  (”2  *  1"3)ltM  *  “T ♦  *• 


J*ki 

hr 


nix2  *  n2n3^3 . ,  nr* 


n2n3X2 


•  “3  M  VK? 


Bov  if  Equations  (13  -  lfc)  are  solved  for  K-^  and  K^, 
*DA  -  *D  -  *1 


hm*L+*B 


<F2r  +  *V 


)  -  Kin,  -J.  -  n3%)  ^ 


-  fag  +  ■  -^L3^)  -  1(1^  -  n3%) 

+  png  +  yaffil-3)  +  1(6^  ♦  i/1  -  n3^3) 

L  V1  *  n32  '  •  J 


vKp  (D6) 


Site 


•  / . 

„ 

/ 

! 

. - 

t 

■  '  ■■] 

/■  -  -  ■  , 

•  • 

In  order  to  eiwplify  Equation  (D 6)  we  will  consider  only  those  terms 
which  are  cubic  or  less  in  or  n^.  Since  Equations  (62  -  83)  indicate 

dx2  ^3  A 

•jjr  and  -g~  are  second  order  in  r.g  or  ny  Equations  (66  -  671)  for  Xg 


and  X,  reduce  to 

t.  -  1,2 

^a/TTT 


1  -  °lB; 

'-.ITT. 


But  both  B  and  J  are  multiplied  by  an  n£  or  n^  wherever  they  appear 
in  Equation  (D6)  ,  and  hence,  we  need  only  consider  second  order  approxi* 


atlons  to  them. 


2  .  _  2 


•  *2  _  *  2  .  *  2  n2+nln3  .  2  2 

.  .  6  -  X,  +  X3  - - -  ^  y  -  ^  +  n3 


-  8“  *  n^ 


(D  10) 


Inserting  Equations  (D7  -  D  10)  in  Equation  (D6)  and  regrouping. 


,  dXj,  dx,  p  dx„  dx  ” 

(F2r  +  i?3r)  "  “  ?  ^)KD  +  (n2  _  dZ[)  +  l(n3  ' 

I*  — 


-  -  +  +  -  “3  >  (Kj,  -  ivKy)  (D  U) 

J 

Revton's  equation  for  motion  in  the  XgXj  plane  may  be  written  in  the 
£orm  ®(*2+  ix^)  «  F2r  +  iFJr  -  img  -  aae  (D  12) 


where  ■  la  mass 

*G  is  the  Coriolis  acceleration.* 


An  expression  fra*  ac  is  given  in  Reference  5. 


i  in ji  I  1 1 Tin l  ii  n  f i |>i ii t  1 1 r w.>>- . ; .y.^. 


? 


If  ve  change  the  Independent  variable  in  Equation  (D  12)  from  tine  to 

t 

die^nsionlesB  arclength  p  »  f  §it,  and  solve  for  the  acceleration, 


V  +  W  /*2  +  SxU'  d  ,  i-  x  .ed  aC4  (r 

- 5 - ■  *  v — ^ - >3  +  —j  (F2r  +  1P3r)  "1T  "  ~T  13 ) 

mu  u  vi 

Substituting  in  Equation  (D  13)  fioa  Equation  (D  11)  and  noting  that 
drag  ter*  in  (D  11)  differs  from  the  ^  term  in  (D  13)  by  a  fourth  order 
term,  v*  can  Integrate  Equation  (D  13)  and  obtain  the  folloving  form 


(x2  +  ^  *  4<*3  +  *30  +  x30}  *20  +  ix30  .  ^0  +  ^30  _  .  ^  ^ 

- j - * - 5 - + - 3 - p  +  SL  '  1SF 


uhere  3^  -  iSj, 


A  |  -  “S’* 


o  1/2 

(1  -  n3  )  dr  dq 


2  "2 


T  .£*-  r 


«d2  ^  ^  u"  2  dr  dq 


1  L 

2  -2 


2  ’  2 


Since  nj  i  (t  - 

idm*  I  -  lXj 


«  2 

ac  u  dr  dq 


■  C.  ! 


-r - i_ 


(1  -  n5  ) 


1  2 

■*■  *  2  a3  +  '  ‘ 


1  +  T!  0-  ~  “^)2  +  •  •  •  » 


«  i  -  ^  52  +  i  (l2  +;S2)  +  .  .  . 


(D  15) 


If  the  assumed  eplcyciic  solution  for  were  substituted  in  the  third 
tern  on  the  right  side  of  Equation  (D  15)  a  simple  calculation  would 
show  that  the  resulting  expression  does  not  contain  either  a  constant 
term  or  terns  in  the  difference  frequencies.  This  term  In  Equation  (D  15) 
will,  therefore,  he  placed  in  our  error  expression.  If  a  linear 

p 

dependence  of  the  lift  coefficient  on  8  is  assumed  and  fourth  powers  of 
K10  's  are  neglected, 

sl  -  S’  J  +  *2  +  vhK  ^  +J*2]  *  d*+  -D*16) 

X/  L  "" 

2  ‘5 

where  -  J^2  -  (JjJ^ 

Inserting  the  general  expression  for  a  damped  epicycle  for  \  in  Equation (E  16) 
end  regrouping, 

^  "  L  ^  *1  '101  *  [\  *  V  *<*  ]  "*  ♦*! 

"  2  -  2  L 


O  17) 


8el2  *  K102  +  ^O2 


ae22  *  +  ho2 


p  q  f  i 

h‘  \  J  |%2  V  *  <B>JL  ^  «*  4* 


L  L 
2  '  2 


HffliWunftHu  iiSilfcr  jiiife 


ii  j) 


f"  2  pi  2 

|<K:'  *  4  -<V  -  Ki'u  j  K‘i  e  4  |_2(Ki  -  Kic 

P  ,]  o 

4  <K2  -Km  )J  K2c  4ltrK?e 


CiKllarly, 


+K1K2  6 


1(2(9.,  -  fU 

<~  A. 


?  Q 

1  NIVv^ 

!  ‘  S  L 


51  .  .li(l 


r  2]  i^2~' 

4  K  4  V  *e2  J  K2  *  ^  +B/ 

p  <i  r  1 

*r  •  ■  S  5  [_v  VEa5^V  Jf<*2  +^J dq 


L  L 

*  7  ‘  S 


(I)  18) 


V  *  V  *  *  % 

Equation*  (0  17  -  0  18)  may  t>*  used  either  for  the  analysis  of  the 
resultn  of  present  swerve  reductions^  or  for  the  construction  of  a  new 
swerve  reduction.  If  the  el gat  unlnovn  epicycle  swerve  reduction  Is 
considered,  Equations  (D  17  -  0  15)  without  the  error  terns  indicate 
that  the  lift  and  Magnus  force  coefficients  as  obtained  fro*  a  given 
■ode  of  oscillation  should  be  plotted  against  the  approprlat  effective 
yaw  squared.* 

It  la  interesting  to  note  that  for  a  rapidly  spinning  model  although 
only  the  swerve  assoc leted  with  the  slow  are,  Kp,  .*  as as ur able,  the 
squared  amplitude  of  the  fast  arm  has  twlcs  as  muclv  wslght  In  the 

p 

calculation  of  6^  as  that  of  the  slow  arm. 


The  mors  modern  LF  sverve  reduction  any  be  applied  to  w.  *>*uur4wnt 
of  non-31rear  forces  when  the  model  Is  rapidly  spinning.  Only 
the  slow  arm  maxes  s  measurable  contribution  and  the  coefficients 
obtained  from  this  reduction  should  be  plotted  against  6  S. 


/ 


,  v.T-0^  8*e-ve  reaction  may  be  constructed  vhtch  vtli  not 

A  special  y^r^e.  sverv,  r 

«  _e  *e<~  *0  ve  rewrite  Ec’iatlons  [  D  17 

-\ect  'ch<*  error  terns.  <-o  «**■*» 


SL  ’  fJL0  +  V  (K102  +K^U2)]^  S 
L  *  |  ‘  2 


^  dr  dq 


(D  19) 


5  £  P  155,  2 

v  S  S  ^  Ki'  *•*  1 


dr  dq  +  E. 


°  h  L 

""5*5. 

[w  ?  2 
JFQ  +  JF02  (K10  *  K20 


t)  n 

in 

'  i  “2 


'A  dr  dq 


(D  20) 


P  4  10  p 

^  1*Kl°  **' 


L  X- 

5*5 


If  estimate*  are  made  for  JF0  ^  JFb2'  **  1 

the  right  aide.  of  Equation*  (D  19  -  D  20)  «  -own  functions  of  P  «A 
^  *«  aubstracted  from  the  measured  J^—2  of  Equation  (D  It,  If 

Z  modified  Equation^  It)  in  .1*  !•  ^  *  looot  tquarcs 

to  the  measured  Xg  individual  round  values  of  ♦  ^2* 


*2+1* 


SW 

5/ s,  individual  round  values  of  Jj^  ♦  ^2 


2  2,  .  *5*  rv24-K2)  may  be  obtained.  When  these 

(Kl0  +  *20  >  *n4  JF0  +  JF62  K10  +  K20  ' 

for  dlfforont  kMlo  of  tt»  »“»  oof  lI«r«Uoo  ««» 

»oh  ™«r  or.  pl««d  win.*  K102  *  *»’ «“  '°™ 

.  thev  dlfftr  from  the  Initial  guaases 

coefficiant#  may  ba  obtained.  Should  they  oirrtr 

th.  proc.se  can  be  Iterated.  This  reiuction  ha,,  been  coded  for  the 
ORDVAC  and  Is  available  for  use. 


► 


APPENDIX  E:  LARGE  YAWING  MOTION  WixH  ALL  30E  LINEAR  MOMEITIo 

In  this  Appendix  we  now  wiU  apply  our  techniques  to  the 

equation  with  non-linearities  in  all  the  moments.  Only  a 
dependence  of  the  moment  coefficients  on  a  will  be  considered. 
«  1  ,.2,'  T 


respectively,  that  equation  may  be 

written 

X"  ♦ 

*  2 

+  Bg  8 

+  |  (62)'  - 

lij  X 

- 

[f 

*  2 

Mg  6 

-  JT  2(&2)' 
“6 

♦  iv(T0 

where  KQ  - 

V 

V 

*  J 

rMA0) 

* 

H2“ 

V 

'V 

+ 

'  H2 

m 

* 

*  V 

-2  * 

+  k2  V 

*  T  1  t 

Ji-  V?  \ 


V '  V  'V*®  ,“o 


_  -2  . 
k2  JM 


Substituting  Equation:  (24-)  for  a  damped  epicycle  in  Equation  (El)  and 

-•'ting  the  coefficient  of  K.e^l  to  zero,  ve  obtain  the  following 

•  i  2 
— .ion .  (Our  basic  assumption  requires  that  K.  ’s  which  should  be 

*•  p 

in  Equation  (E2)  be  replaced  by  midrange  value  KiQ  .) 


-j 

S  +  M26el2)  +  i7(T0  +  T2  *el2) 

+  4  f^O2  +  K202)  °1  +  +  *20^  +  ^2^ 

•  QU  +  P  o 

■  ^1^2  +  5  )  ('2al  *10  +  2a2  ^20  ^ 

*<%2*~ 4—  )K202  ^  "  t4>j  - 
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A  similar  equation  can  be  derived  from  the  coefficient  of  K^e  .  If 
squares  of  Ao  are  aaitted  in  the  equations  for  (  -  +  10| )  and 

(-  ♦  10?),  the  familiar  relations  for  the  linear  case  may  be  obtained; 


"o  -  *[  •  <  *  «w» 


=0  -  °1  +  «fe 


‘  <+*2  ' 


81nce  in  the  linear  case  the  equation:  for  V  is  unaffected  by  damping 

and  that  for  Mq  has  a  very  small  term  from  the  damping,  we  would  expect 

the  addition  of  damping  to  have  little  effect  on  Equations  (72)  and  (93) • 

If  we  make  the  simplifying  assumption  that  this  Is  the  case,  the  real 

1.0 

part  of  Equation  (E2)  and  its  equivalent  for  e  2-wuld  yield  very 

slightly  modified  versions  of  Equation  (92)  and  (95)  and  only  the  Imaginary 

# 

parts  of  these  equations  need  be  considered. 

(EV 

““X  -  Vl  -  *  «i  ‘  *<*0  *  4  ».l‘>  -  4  [K10S  -  «202*1  *  **,% 

The  precise  effect  of  the  damping  is  the  addition  of  a^(H  -  a^)  to 
lq.  (91)  *nd  (^(H  -  0^)  to  Eq.  (92)-  (See  Eqs.  109-110). 


-fcti 


»  •  p  t  o  O  p 

"  *  ^*20  +  ^alK10  +  “2^20  ^  +  ^  2  ^20  "  0 

*  r  '  'i  (EI5) 

Vs  -  «*>  *  »(I0  *  T*  »./)  -  B*  [(iCj/  .  K,J)('Z  *  K^J 

-  \i  <4  -  4>icio2 ♦ *  °y^o2»  *  ^K—Ko  - 0 

o 

The  appearance  of  6g^  in  Equations  (E4)  and  (E5)  could  have  been  predicted 

from  the  fact  that  vT  is  the  imaginary  part  of  the  coefficient  of  X 

while  M  la  the  real  part  of  the  same  coefficient.  Thus  for  a  Magnus  force 

o 

which  Is  a  polynomial  function  of  6  ,  the  relations  in  Table  1  may  be  used 
in  the  generalized  forms  of  Equations  (E4)  and  (E5)  in  the  same  way  that 
they  may  be  used  for  a  polynomial  static  moment  in  generalized  forms  of 
Eqa.  (92)  and  (93). 

T-  the  small  geometric  non-linearity  terms  are  omitted*  linear  damping 
moment  i  are  assumed*  and  the  effect  of  the  small  non-linear  lift  terms  is 
neglected*  Equations  (EU)  and  (E5)  assume  the  simple  forms 


Va  -  »(*Q  *  <  *J) 

20o  -  V 


An  examination  of  Equations  (92)  and  (95)  shows  that  when  the  ,,Jomotric  non¬ 
linearity  terms  are  neglected*  the  frequencies  depend  on  the  amplitudes 

2 

through  their  corresponding  effective  squared  yaws,  !>el  .  Thus*  ve  have 

the  important  result  that  under  these  assumptions  the  damping  of  a  given 

mods  Is  related  to  the  amplitude  of  the  total  motion  by  means  of  Its 

effective  square  yaw.  In  other  words,  for  fixed  aerotvnemic  characteristics 

and  spin,  the  damping  exponents,  a,*  are  functions  of  tnei*  *  2,s  alone. 

z  ox 


4,  ' '  ■£  *  yy  ■ 


\ 


1 1* 


If  Equation  (E4)  is  now  subtracted  from  Equation  (E5)>  the  result 
•  solved  for  H0  and  simplified  by  use  of  Equation  (9*+) 

*  -  0.^1  q2.  -  *  ,K10  '  K2o\ 

=o  *  A  w  • ’*« 


*  (Oj-OjK^q  -KgQ  )  2  2^ 

M-  —  y  ~779 -  '  \2  (K10  ^0  ’ 

*  (0X  -  02  r  6 


vhere*  ®rsnge  *  «i  +  «2  + 


(0X  -  ~  (02  -  ^l)g2KgQ- 

2(0X  -  _  . 


Eliminating  Hq  from  Equations  (E4  -  E5)  and  simplifying  by  use  of  Equation  (9,), 
and  ignoring  all  K,  terms,  we  obtain  the  following  equation. 


T0  ♦  T2  “«2  *  W 


.  .\o%2  '  >WV \ 


where  T 


.  <«,  -  ^Ho  -  *co2>  ,  *  *«L*g) 

“* - h 

.  [1  -  (Ix^i ' 

01  +  02  [  h  -  *2 


0:0l(°h-°b)(Kio2^  ;  "aX> 

(0X+  0^  .(0!  - 

The  quite  formidable  expressions  for  Hrange  and  TranGe  differ 
appreciably  from  the  simple  relations  for  Hq  and  fQ  of  Equations  (E3) 
only  when  the  yaw  is  qpite  large. 

•  The  subscript  "range"  denote  the  quantity  which**  range  reduction 
TOUld^et  if  only'the  geometric  non-linearities  were  considered. 


defined  by  the  relations 


The  range  values  of  J„  -  J„*  8111  Jij  83:6 
following  Equation  (20) 


-  _  V  C  T  -  T 

“  Ki  "L«  range 

range  0 


vJH  -  JMA)range  "  V  grange  “  %  '  Vj 
[\  '  Traa?C] 


(E  10) 


(E  11) 


If  these  relations  are  used  in  Equations  (E8  -  E9)  ®0’  ^2’  ‘O'  2 

and  K,  are  replaced  by  their  definitions  from  Equations  (El)  and  (9), 
and  the  J^s  converted  to  K^s  the  following  two  equations  for  (Kg  - 

and  1L,  may  be  obtained.  (E  12) 

*range  «  ?  2 

*  /^lIC20  ‘^0_\ 

•••  (kh  -  K«*,r„8.  •  \  •  X.  *  V  *[ -h 

-  «r  2  (I)(7-Z  ,(K1«*  •  K2°2) 

A  0x  -  02 


+  \2  a2  ^2  +  R5  V  *  %  % 

*  2  *  **1  '  nfO  jg - )  (E  13) 

Sp  “  V  +  *T.2  8e  +  \2(B)(  m'2_  <$} 

Trange  OB  °  ”l  r2 


♦  “1 V  ’  V  ‘  b5  v  * b'  \ 

vhere  "  ag)^^io  “  ^20  ^^1  ”  ^ 

^  .  t,2  -  0X,2  -  -  <)*202]  <*i  *2>  1 
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o  IT*  ?  .»  o  ,»  .  -1 

a3  *  *  k2  j^l  ^0  "  K10  ^1  *  ^ 


%  *  (|)k22  <**“  +  O 


bl  *  I  "  ^2^*10  ‘  ^  ^1  "  ^2  ^ 

b2 . .  tl2  j^;2 . 42)Kl02 .  t^;2 .  <2)kJ  w;2 .  4V1 

>3  -  -  ^  <"io2  ■  *J  &<■*?  -  42*'1 

\  “  ■  &kl2  (K102  *1  +  K202  *Wl  +  4)"1 

The  terms  In  and  b^  in  Equations  (E  12  -  E  13)  are  usually-  smell 
and  can  be  neglected.  (For  those  cases  where  these  terms  are  important 
the  terms  must  be  calculated*  and  subtracted  from  the  range  values  of 
kh-k ma“4V)  Equation  X  2  12)  predicts  a  dependence  or  (hg  -  K^) 

ra 

#  * 

on  both  L  j  and  Kj,  g  while  Equation  (E  13)  predicts  a  dependence  of 

^  *8  ra 

on  both  K^g  and  g. 


Only  three  ro’tnds  of  the  126  round  e  of  the  Army  -  Havy  Spinner  Rocket 
Program  required  this  calculation.  For  one  of  tteee,  however,  the 
correction  to  Kg  -  K^  due  to  a,  K^.  2  was  about  3'*^, of  the  uncbrrected 
value.  *  o 
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